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Abstract 



Within the universality class of ferromagnetic vector models with 0{n) 
symmetry and purely dissipative dynamics, we study the non-equilibrium 
critical relaxation from a magnetized initial state. Transverse correlation and 
response functions are exactly computed for Gaussian fluctuations and in the 
limit of infinite number n of components of the order parameter. We find 
that the fluctuation-dissipation ratios (FDRs) for longitudinal and transverse 
modes differ already at the Gaussian level. In these two exactly solvable cases 
we completely describe the crossover from the short-time to the long-time 
behavior, corresponding to a disordered and a magnetized initial condition, 
respectively. The effects of non-Gaussian fluctuations on longitudinal and 
transverse quantities are calculated in the first order in the e-expansion (e = 
4 — d) and reliable three-dimensional estimates of the two FDRs are obtained. 
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I. INTRODUCTION 



Consider a ferromagnetic system with vector order parameter $j(x, t), % = l,...,n, 
a critical point at temperature T c and at vanishing external magnetic field h c = 0, and 
prepare it in some initial configuration (which might correspond to an equilibrium state at 
a given temperature To and magnetic field ho). At time t = bring the system in contact 
with a thermal bath of temperature T& ^ T in the presence of a magnetic field hi,. The 
ensuing relaxation process is characterized by some equilibration time t cq so that for t ^> t cq 
equilibrium is attained, whereas for < t <C t cq , the evolution depends on the specific initial 
condition. Upon approaching the critical point T b = T c , hj, = h c = the equilibration 
time diverges and therefore equilibrium is never achieved. One-time quantities (like mean- 
energy, magnetization, etc.) typically display a slow relaxation (power-law) towards their 
equilibrium values, even in the extreme situation in which t cq is infinite. Consequently, 
the natural objects to monitor the non-equilibrium evolution of the system are two-time 
quantities, like the time-dependent correlation function of some local observable O which 
is given by C®(t,s) = (0(x, t)O(0, s)), ((•••) stands for the average over the stochastic 
dynamics) and the linear response R® (t, s) to an external field h conjugate to O. R®(t, s) 
is defined by R®(t,s) = 5(0(x.,t)) / Sh o (s)\ ho=0 , where h a is applied at time s > at the 
point x = 0. For t > s ^> t cq correlation and response functions take their equilibrium forms 
- depending only on t — s — and they are related by the fluctuation- dissipation theorem 
(FDT) 

T b R°(t > t eq , s > t cq ) = d s C°(t > t cq , s > t cq ) . (1) 

This is no longer true in the so-called aging regime: correlation and response functions have 
a non-trivial (homogeneous) dependence on both t and s, even for very large times. In 
particular their decays as functions of t become slower upon increasing s, the age of the 
system. This is why this evolution is referred to as aging. 

The FDT suggests the introduction of the so-called fluctuation- dissipation ratio 
(FDR) [1,2]: 

X°(t s) - TbR ^ t ' s ) (2) 

which can be used in order to understand the exten to which a system is evolving out of 
equilibrium. The asymptotic value of the FDR 

X™ = lim \imX° =0 (t,s) (3) 

>oo t— >oo 

is a useful quantity in the description of systems with slow dynamics, since Xq = 1 whenever 
the aging evolution is interrupted and the system crosses over (for t > s ^> t cq ) to equilibrium 
dynamics. Conversely, X@ ^ 1 is a signal of an asymptotic non-equilibrium dynamics. 

Within the field-theoretical approach to critical dynamics it is more convenient to focus 
on the behavior of observables in momentum space. Accordingly, hereafter we mainly con- 
sider the momentum-dependent response R® (t, s) and correlation C®(t, s) functions, defined 
as the Fourier transforms of R®(t,s) and C®(t, s), respectively. In momentum space it is 
natural to introduce a quantity with the same role as X°(t,s) [3]: 

T b B°(t,s) 

ql ' j d s C°(t,s)- {A) 
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It has been argued that the long-time limit 



voo 
A 



= lim lim Xz=q 



(t,s) 



(5) 



has exactly the same value as Xq [3] (see also Ref. [4]). 

Aging was known to occur in disordered and complex systems (see, e.g., Ref. [5]) and 
only recently attention has been focused on simpler critical systems such as ferromagnets, 
whose universal behaviors can be rather easily studied by using powerful theoretical methods 
and which might provide insight into more general cases [6] (see Ref. [7] for a pedagogical 
introduction). In particular it was argued that X°° is a universal amplitude ratio for quenches 
to T c [8]. In this respect several dynamic universality classes were investigated and the 
scaling forms the of response and correlation functions, together with the associated FDRs 
were calculated by means of exact, approximate, or numerical methods [9-30] (see Ref. [6] 
for a review). 

In this paper we consider the purely dissipative dynamics (Model A, according to the 
classification of Ref. [31]) of the 0(n) symmetric model relaxing from a magnetized state 
(in a sense that will be specified below). Previous works [32,33,30] mainly focused on the 
case n — 1 - - corresponding to the Ising universality class - - because the dynamics of 
actual systems with vector order parameter and 0(n > 1) symmetry typically possesses 
conservation laws which are not accounted for by Model A. For example, an isotropic fer- 
romagnet (antiferromagnet) is correctly described by the so-called Model J (G) universal- 
ity class [31]. Nonetheless, we consider this unphysical model because it is probably the 
simplest example in which the dynamics is affected by two different fluctuation modes. 
In fact, consider the relaxation from an initial state with a non- vanishing magnetization 
M . In the absence of external fields the magnetization will evolve as M(t) in the direc- 
tion Mo = Mo/Mo. The time evolution of the component of the order parameter parallel 
to M(t) is clearly different from the evolution of the n — 1 components perpendicular to 
it, which are characterized by a residual 0(n — 1) symmetry. In what follows the mode 
er(x, t) = $(x, t) ■ M — M(t) parallel to M is referred to as longitudinal, whereas the n — 1 
perpendicular modes 7r(x, t) = $(x, t) — M [M • $(x, t)] are referred to as transverse. They 
are characterized by different responses and correlations, with FDRs: 



In spite of its unphysical nature, the purely dissipative dynamics of the 0(n) model has 
a rich behavior which is easily accessible to numerical simulations and possibly captures 
some of the features of more realistic dynamics. The model has already been studied by 
Fedorenko and Trimper [24] in a dimensional expansion close to the lower critical dimension 
— an approach that is complementary to the one adopted here. 

The paper is organized as follows. In Sec. II we introduce the model and we discuss 
the general scaling forms for the correlation and response functions of longitudinal and 
transverse modes after a quench from a magnetized state. In Sec. Ill we solve the model 
within the Gaussian approximation. In Sec. IV we derive the exact transverse response and 
correlation functions in the limit of infinite number of components of the order parameter. 
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In these two exactly solvable cases we describe completely the crossover from the short-time 
to the long-time behavior, corresponding to a disordered and a magnetized initial condition, 
respectively. In Sec. V we account for non-Gaussian fluctuations up to first-order in the 
e-expansion (e = 4 — d) and from the resulting expressions of the response and correlation 
functions we calculate the associated FDR. Finally in Sec. VI we summarize and discuss our 
results. 



II. THE MODEL 



The purely dissipative dynamics of a n-component field $ = (<E»i, . . . , $ n ) can be specified 
in terms of the stochastic Langevin equation 



5H\§\ 

d t ^t) = -n-PP- + &(x,t), 



(7) 



where Q is the kinetic coefficient, £j(x, t) a set of zero-mean stochastic Gaussian noises with 

(&(x,^(x',0> = 2V}<f(x-x')<S(t-t'), (8) 

and is the static Hamiltonian. Near the critical point, 7i[<&\ may be assumed of the 

0(n)-symmetric Landau-Ginzburg form 



nm = f 



d d x 



(V$) 2 + ^$ 2 + |($ 2 ) 2 



(9) 



where is a parameter that has to be tuned to a critical value ro )C in order to approach the 
critical temperature T = T c (r 0)C = 0, within the analytical approach discussed below), and 
go > is the bare coupling constant of the theory. 

Correlation and response functions of a field which evolves according to the Langevin 
equation (7) can be obtained by means of the field-theoretical action [34,35] 



S[$, $] 



dt / d d x 



5$ 



(10) 



where $(x, t) is an auxiliary field, conjugate to the external field h in such a way that 
W[$,h] = H[$] -J d d xh-$. 

The effect of a macroscopic initial condition $o( x ) = ^( x , t = 0) may be accounted for 
by averaging over the initial configuration with a weight e~ H °^°\ where [36,37] 



H [$o] = y"d^^[$o(x)-M ] ; 



(11) 



specifies an initial state with Gaussian distribution of the order parameter fluctuations, 
short-range correlations proportional to r^ 1 , and spatially constant average order parameter 
M = ($ (x)). In the experimental protocols we have in mind the initial state with M ^ 
can be obtained by equilibrating the system for t < in the presence of an external field h 
and at generic temperature T . The resulting state can be described by Eq. (9) with external 
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field ho, leading, far enough from the critical point, to a Gaussian distribution of the order 
parameter such as Eq. (11) with r ~ goMfi where M = (6|h |/(7o) 1//3 - However, within the 
renormalization-group (RG) approach to the problem it has been shown [36] that Tq" 1 is an 
irrelevant variable in the sense that Tq 1 affects only the correction to the leading long-time 
scaling behavior we are interested in. In view of that we fix it to the value r _1 = from the 
very beginning of the calculation. It is important to note that the initial state described by 
by the Hamiltonian Eq. (11) does not correspond to a low-temperature one in zero magnetic 
field, unless n = 1 [30]. In fact, for n > 1 the transverse modes are critical for all T < T c 
(Goldstone theorem, see, e.g., [34]) and the proper H takes the form of Eq. (11) only for 
the longitudinal mode, whereas the most relevant term (in the RG sense) for the transverse 
ones would be J d d x (V7r ) 2 (where 7To(x, t) = 7r(x, t — 0)). 

In order to account for a non-vanishing mean value of the order parameter ($(x, £)) = 
M(£) (we assume that M(£) stays homogeneous in space after the quench) during the time 
evolution, it is convenient [33] to write the action (10) in terms of the transverse and longi- 
tudinal fluctuations (ir and a, respectively) around M(t), i.e. 



$(x, t) = (<t(x, t) + M(t), tt(x, t)) , $(x, t) = (cr(x, t), tt(x, t)) . 



(12) 



With this parameterization (cr(x, £)) = 0. The problem we shall consider below is there- 
fore the dynamics of fluctuations of the fields (cr(x, t), 7r(x, t)) in the time-dependent "back- 
ground" provided by M(t). For later convenience we also introduce the scaled magnetization 
m(t): 

M 2 



m =9o—, 

so that a perturbative expansion in go leads to a finite value for mo 
resulting action in terms of a, a, tv, tt may be written as 



d d x(C + A + C 2 ) , 



(13) 

m(t = 0). The 
(14) 



where 



C = -atta + a(d t + tt [-V 2 + r + m 2 (t)] a) 

m 2 {t) 



d = n 



-V" + r + 



V2m(t) 



7T, 



—era H — — or H — (T7r7r I 

O O J 



. 9o ~ 3 , 9o ~ 2 . 90 ~ 2 , 90 2 ~ 
H o~ o~ H aair H irrnr H a tz7t 

6 6 6 6 



a a + n 




m(t) = ah eS (t) 



(15) 

(16) 
(17) 



We split up the action S so that Co is the Gaussian part, C\ contains the interaction 
vertices and C2 gives the coupling to the effective magnetic field h e g(t), due to a nonzero 
M(t) oc m(t), acting on cr(q = 0,t). The presence of m(t) 7^ breaks the original 0(n) 
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symmetry of S in the fields $ and $ (see Eq. (10)), leaving only a reduced 0(n — l) symmetry 
in the transverse field iz and tt (to all orders in perturbation theory). 

Following standard methods [34,35] the response and correlation functions may be ob- 
tained by a perturbative expansion of the functional weight e -( 5 '[*>*]+- f/ o[ < i > o]) j n terms of the 
coupling constant g . The propagators (Gaussian two-point functions of the fields a and 7r 
in momentum space) are [36,37] 



(a(q, t)a(q!, f )>o = (2n) d 5( q + q')R^(t, f ) , 
(a(q, t)a(q', f)>o = (27r) d 5(q + q')C£(f , f ) , 
(^(^0^(^,0)0 = (27r)%5(q + q')^(t,t'), 
(^(^0^(^,0)0 = (27r)%5(q + q')Q(t,t'), 



where (g = |q| 



i^(M') = 6(t - t') exp <^ -n 



(q 2 + r )(t-t') + f dt"m 2 (t") 
Jf 



R£(t, = 0(t - exp (g 2 + r )(t - *') + J* dt"m 2 (t")/3 j , 

/•OO 

CJ w (t, f') =2Sl dt"R^ n (t, t")R a ^{t', t"). 
Jo 



(18) 
(19) 
(20) 
(21) 



(22) 
(23) 
(24) 



The perturbative expansion can be as usual organized in terms of Feynman diagrams. The 
response propagators are given by Eq. (22) for the longitudinal mode and by Eq. (23) for 
the transverse one, and they are represented as directed solid and dashed lines, respectively, 
where the arrow points towards t, the larger of the two times. The correlation propagators 
are given, instead, by Eq. (24) and are represented as solid and dashed lines for longitudinal 
and transverse modes, respectively. The interaction vertices are contained in L\. Four of 
them are standard time-independent quartic vertices oc flg whereas the remaining three are 
time-dependent cubic ones oc Q,^Jg^m(t), due to a non-zero magnetization m(t). Propagators 
and vertices are shown in Fig. 1. The contribution of h e e(t) cancels for a suitable choice of 
m(t). The evolution equation for the magnetization m(t) can be obtained by solving the 
equation of motion (5S/5a(x, t)) \a=o = 0: 



d t + n [ r + -m 2 (t) 



O 3//2 

m(t) + Q^-m(t)(a 2 + tt 2 /3) + -^(a 3 + aiz 2 ) = . (25) 
2 6y2 



To the lowest order [hereafter (dq) = d d q/(2it) d and VL = 1] 

( ( x 2 (x,t)) = J(dq)C°(t,t) + O(g ), ( ( x 3 (x,t))=0( v ^), 
(tt 2 (x, t)) = (n - 1) J (dq)C^(t, t) + O(g ) , (<t7t 2 (x, t)) = 0{^) , 
and therefore up to one loop and at the critical point r = 0, Eq. (25) becomes 



= d t m(t) + 



+ 



(dq) 



c°Jt,t) + \±cz(t,t) 



+ 0(g 2 ) 



(26) 
(27) 

(28) 



6 



R^t,s) C°(t,s) 
> 



--> 



50/6 50/6 50/6 50/6 



-><- 

FIG. 1. Elements of the diagrammatic representation of the perturbation theory. 

A. Scaling forms 

When a critical system is quenched from the high temperature phase to the critical point 
response and correlation functions display the scaling forms [36,6] 

i? q=0 (t, s) =A R (t- s) a (t/s) e F R (s/t) , (29) 
C q=0 (t, s) = A c s(t - s) a (t/s) 9 F c (s/t) , (30) 

where a = (2 — r\ — z)/z, z is the dynamical critical exponent, 77 the anomalous dimension 
of the field, and 9 is a genuine non-equilibrium exponent [36] . Ar and Ac are non-universal 
amplitudes which are fixed by the condition -Fr,c(0) = 1. With this normalization Frc are 
universal and X°° turns out to be a universal amplitude ratio X°° = Ar/[Ac(1 — 0)\ [8,6]. 
Although we shall focus mainly on the case q = 0, the generalization of these scaling 
forms to non-vanishing q amounts to the introduction of an additional scaling variable 
y = Aci$lq z (t — s). An is a dimensional non- universal constant which can be fixed according 
to some specified condition. This scaling behavior gets modified when the critical system 
evolves from a state with M 7^ 0. The magnetization scales according to [36] 

m(t) = A m m t a+e F M (B m m t K ) (31) 

where k — 6 + [35/ (vz) — 1 and standard notation for critical exponents has been used. 
The non-universal amplitudes A m and B rn can be determined, e.g., by imposing JF M (0) = 1 
and = — 1 [30]. From Eq. (31) one sees that the effect of a non- vanishing initial 

magnetization itlq is the introduction of an additional macroscopic time scale r m into the 
problem, i.e., r m = (_B m m ) _1 / K and of an additional associated scaling variable u = t/r m 
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in the scaling forms. This scaling variable, naturally appears in the scaling forms for the 
response [Eq. (29)] and connected correlation functions [Eq. (30)]: 

i? q =o(*> s ) = A * (* - s ) a (*/s)"^R(s/*> Sm^o*") , (32) 
Cj=o(*, s) = A c s(t - s) a (t/s) e r c (s/t, B m m t K ) , (33) 

where no new non-universal amplitudes have been introduced. The index % = a, tz indicate 
the longitudinal and transverse fluctuation modes, respectively. The resulting functions J-r 
and Tq are universal. For m = the original 0(n) symmetry of the model is restored and 
therefore transverse and longitudinal response (and correlations) become identical. 

According to Eqs. (32) and (33) a non-vanishing mean value of the initial magnetization 
m 7^ affects the scaling properties of the response and correlation function as soon as 
B m mot K ~ 1 (i.e., t ~ r m ) and in particular this happens in the long-time limit we are 
interested in, characterized by t ^> s ^> r m . This formally corresponds to the case m — > oo, 
as opposed to the case previously considered, m = 0. In this limit one expects the scaling 
forms (32) and (33) to turn into: 

R^o(t, s) =A%(t- sy(t/s) e °F°(s/t) , (34) 

C° =0 (t, s) = A° c s(t - sY(t/s) e '°F° c (s/t) , (35) 

Rl =0 (t, s) = Al(t- s) a (t/s) e -F-(s/t) , (36) 

C- =0 (t, s) = AZ s(t - s) a {t/s) e '-F-{s/t) . (37) 

The non-universal constants A RC are fixed by requiring F RC (0) = 1, where F RC (x) are 

universal functions related to the large-f behavior of J- RC (x : v). The exponents Q, L and Q\ 
are different from 0. In particular it is not obvious, a priori, whether 9{ and 9[ should be 
expected to be the same, if they are novel exponents — as 9 — or just combinations of 
known ones. 

In Ref. [30] we showed, for the Ising model (i.e., the 0(n) model with n = 1), that 

: = = (38) 

The RG proof of Ref. [30] made no use of the scalar nature of the order parameter and holds 
for any purely dissipative model (i.e., with Model A dynamics). 
For the transverse modes, it has been argued in Ref. [24] that 

^ = ^ = ~- ( 39 ) 

A fundamental difference between the quench from a disordered and a magnetized state 
emerges: In the latter case the aging properties and the non-equilibrium behavior of both 
longitudinal and transverse modes are completely described in terms of equilibrium expo- 
nents. 

As a consequence of the scaling forms Eqs. (34), (35), (36), and (37) the two limiting 
FDR are 

X°°(m ^ 0) = - f R - . , A°°(m + 0) = - „ . , (40) 

which turn out to be universal amplitude ratios [30]. 
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III. GAUSSIAN APPROXIMATION 



The starting point for the analysis of the dynamics resulting from the action S [see 
Eqs. (10) and (14)] is the Gaussian approximation in which one neglects all the terms which 
are of degree higher than 2 in the fields 7r and a, contained in L\. The result of such 
an approximation is exact if the space dimensionality is greater than the upper critical 
dimension d u = 4 [34,6]. 

Within this approximation, the equation of motion (28) at the critical point r = 
becomes 

d t m G (t) = - l -m%{t) , (41) 
(G stands for Gaussian) whose solution is 

^ = 3r^ = < a ^ 1 + dy~'- (42) 

This expression agrees with the scaling behavior (31), given that 6 = a = and k — \ within 
the Gaussian model. The non-universal amplitudes appearing in Eq. (31) are A m = 1 and 
B m = a/2/3, leading to r m = (_B m m )" 1//K = 3/(2mg) and m G (t) = m (l + t/r m ) -1 / 2 . For 
* > r m , m G (t) ~ (2t/3)" 1/2 and a nonzero m affects only the corrections to this leading 
long-time behavior. In this sense 1 is irrelevant for long times. 

Gaussian response and correlation functions (see Eqs. (22) and (24)) do not depend on 
n and therfore the longitudinal ones are the same as those of the Ising model studied in 
Ref. [30], which the reader is referred to for details. We recall here that the Gaussian FDR 
turns out to be = 4/5 [30]. For the transverse critical response and correlation functions, 
Eqs. (23) and (24) give (t > s) 

1/2 

-Q 2 (t-s) (43) 



i 

-q 2 (t+s) 



K*fl(t,s) = e-i 2 (*-*)-/>' m G('')/3 = ( S + Tm \ 

\ t + T m J 

C^(t,s) = 2 dt'Rq ,0 (t,t')Rq'°(s,t') = 2-—— — — -r— / dt'(t' + r m )e 2 ^. (44) 

JO IV' ' 'm)\ s i im)\ 1 JO 

In particular, for q = one finds 



1/2 



^Htt^J • (45) 

Comparing these results (for r m = 0), with the scaling forms (36) and (37) we can identify 
z = 2, a = 0, 9k = — 1/2 in agreement with standard mean- field exponents {v — ft — 1/2 
and r\ — 0) and determine = 1, A £ = 1, F][(x) = 1, and Fq(x) = 1. Accordingly, the 
asymptotic FDR is given by (see Eq. (40)) 

A% 2 

x " = Az(i-e 1t ) = 3' (47) 
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(a) ^ (b) 

FIG. 2. Transverse FDR X* within the Gaussian approximation (solid line) and in the exactly 
solvable limit n — ► oo for d = 3 (dash-dotted line) and d — > 2 (dotted line in the right panel). Left: 
X£ =Q as a function of s/r m . Right: for r m = as a function of g z s. 



which differs from the corresponding quantity for the longitudinal mode X%° = 4/5. To 
our knowledge this is the first case in which the asymptotic FDR of two local quantities 
differ already within the Gaussian approximation. In fact in Ref. [20] we showed that, 
for a quench from the disordered phase to the critical point, the Gaussian Xq is always 
1/2, independently of the specific quantity O it refers to. As a consequence an effective 
temperature T eff = T c /X°° is not well-defined even within the Gaussian approximation for 
the quench from a magnetized state [38,39,20]. 

Let us discuss in more detail the effect of a non-zero magnetization on the FDR 

d,CZ'°(t, s) 

From Eqs. (45) and (46) one finds, as a function oiu— s/r n 



X?(t, s) = q y ! . (48) 



(49) 



which is actually independent of t (a typical property of the Gaussian approximation) and 
increases monotonically from the value 1/2 for u — 0, to the value 2/3 for u — > oo. Ac- 
cordingly, whenever m is non zero, the transverse FDR approaches 2/3 in the long-time 
limit (and its value is indeed independent of the actual ^ 0) whereas the FDR is equal 
to 1/2 for m = 0, i.e., for quenches from a disordered initial state. The plot of X* =0 (s,t) 
as a function of u — s/r m is reported in Fig. 2(a). In the case r m = oo (corresponding to 
m = 0) it was found that the fluctuation mode responsible for aging is the homogeneous 
one q = 0, given that for q ^ and long times X° = X* = 1 [3,6]. The same conclusion 
was drawn for the (longitudinal) fluctuations of the Ising model with r m = (corresponding 
to m = oo) [30]. Now for the transverse modes one finds, from Eqs. (43) and (44), as a 
function of y = q 2 s, 
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which monotonically increases from the value 2/3 at y — to 1 for y — > oo [see Fig. 2(b)], 
as observed in all the other cases studied so far. 



IV. DYNAMICS OF THE TRANSVERSE MODES IN THE LIMIT TV -► oo 

Critical transverse response and correlation functions can be analytically calculated in 
the limit of large number of components of the order parameter. In fact, assuming that 
g = g/n (with g independent of n), the action (14) for large n may be written as (Q — 1) 



S(n,ir;m) = J d xdt |— irn + ir 

where m(t) satisfies the equation of motion 



^+(-V 2 + r + ^ 



7r + f-7f7r7r 2 l+O(n ), 
on J 

(51) 



U + r + ^ + |^ + 0(n-i)) m(t) = . (52) 
\ don / 

It is evident that only transverse modes contribute to the leading term of the action. On the 
one hand, this makes the calculation easier but on the other the behavior of the longitudinal 
mode cannot be determined without including the subleading 0(n°) terms in the action. 

To proceed further we self-consistently factorize the fourth order interaction term ac- 
cording to [34,42] 

-TT7T7T 2 -> C£ =0 (t, t)TZ1T . (53) 

n 

The model resulting from this factorization is Gaussian with a time-dependent temperature: 

r(f) = r + f C^(M) + ^ ■ (54) 
The response and correlation functions are (r^ 1 = 0) 

J R-(t, S ) = exp|-^ < dt'[g 2 + f(t / )]} , (55) 

C%(t, s) = 2 T dt'Rl(t, t')Rl(s, t') , (56) 
Jo 

with the self-consistency condition [r^ = r(t — > oo)] 

r(t) = ^ + | j\dq)[C-(t, t) - CJ(oo, oo)] + ^ , (57) 

where (oo, oo) = l/(q 2 + r^) is the equilibrium correlation function and A is a large- 
momentum cut-off. The equation of motion (52) becomes 

d t m(t) + r(t)m(t) = , (58) 
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which is solved by 



m(t) = m exp 



dt'f(t') 



(59) 



It is possible to cast the previous equations in a very convenient form by introducing the 
function 

h(2A 2 t) = \ exp 



dt'f(t') 



JO 



(60) 



(h corresponds to the function g usually introduced in the literature on the dynamics of the 
spherical model [8,26]) in terms of which 



m(t) = m [2/i(2A 2 t)]- 1 / 2 , 



R"(t,s) =e" q2( *- s) 



h{2A 2 s) 
h(2AH) 



1/2 



(61) 
(62) 



and the expression for readily follows. Accordingly, once we know the function h(x) for 
a generic initial magnetization, the dynamics of transverse modes is completely determined. 
The details of the computation of h(x) are reported in Appendix A. The final result for the 
critical case (f m = 0),2<rf<4(e = 4 — d), and in the scaling regime t ^> (2A 2 ) -1 is 



:t n 1 sin(7r(rf-2)/2) 2 /2 
u* 7rr((d-2)/2) 1 1 



h(2AH > 1) 



1 + 



rrint 



(63) 



where u* = g*N d /(6A t ) is the fixed-point value of the four-point coupling (see Appendix A) 
and corrections to the leading scaling behavior have been neglected. 
According to Eq. (61), the magnetization is given by 



m(t) = m Q 



u 



7rr((rf-2)/2) 1 1/2 
2sin(7r(d- 2)/2) 



(2A 2 t) 



1 + 



3(d-2 



-m\t 



1-1/2 



(64) 



which can be cast in the scaling form Eq. (31) with (we recall that a = 0, /35/(vz) 
(d + 2)/4 [34], 9 = e/4 [36], and therefore « = 1/2) 



it 



7rr((d-2)/2) 
2sin(7r(rf-2)/2) 



1/2 



(2A 2 Y / \ B r , 



3(d-2) 



1/2 



and 



^m(^) = (1+^ 2 )- 1/2 . 
The time scale associated to the initial magnetization is therefore 



(65) 
(66) 



r m = (B m m )~ 1 ^ = -(d-2)m 2 



The response function follows from Eqs. (62) and (63): 



s 



e/4 



1 + t/T„ 



1/2 



(67) 



(68) 
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which displays (for q = 0) the expected scaling behaviors (29), (32), and (36), for r m = oo, 
generic r m , and r m = 0, respectively. [We recall that 6^ — 9'^ — —f3/(vz) = — 1/2 + e/4.] In 
particular, comparing for r m = to Eq. (36), one finds = 1, Fr(x) = I. The correlation 
function Cq (t, s) can be easily computed from Eq. (68) via Eq. (56). In particular, for q = 
one finds 



°q=0 



, \ 4 

(t,s) = -s 



(»"' 


1 + s/r m 


1/2 - 




_ 1 + t/r m _ 





1 



(69) 



d-21 + s/t % 

which displays also the expected scaling behaviors (30), (33), and (37), for r m = oo, generic 
r m , and r m = 0, respectively. In particular Eq. (36) holds with A 1 ^ — 4/d and Fq(x) = 1. 
Equations (68) and (69) allow the computation of the FDR: 



d+2 



1 + 



d 2 -4(l + u) 5 



(70) 



where u = s/r m . For m = (i.e., s C r m oc m o 2 ) we recover the well-known result 
X°° = 1 — 2/d [8,20], whereas for m ^ the long-time asymptotic value for s ^> r m is given 
by 

^(n = oo) = -^-. (71) 



d + 2 



Close to d = 2 this expression becomes 



X?(n = oo) = l+ 1 -e + ()(?) 



(72) 



where e 
therein): 



d — 2, in agreement with the limit n — > oo of the result of Ref. [24] (see Eq. (18) 



X°° = - + 



In- 1 



e + 0(e 2 ). 



(73) 



2 8n-2 

In passing we note that Eq. (71) coincides with the FDR derived in Ref. [26] for local spin 
observables in the spherical model (see Eq. (8.12) therein). 

A plot of Eq. (71) is provided as a dashed line in Fig. 6 (right panel) together with the 
result of the expansion (73) for n = 3 and n = oo (denoted by n = 3[e] and n = oo[e], 
respectively). As in the Gaussian case, A q T =0 (-u) can be used to describe the crossover 
from magnetized to non-magnetized initial state and the scaling function resembles closely 
Eq. (49), — derived for the Gaussian model and corresponding to the case d > 4 — the only 
difference being the values at u — and for u — > oo and the actual expression for 
function of m . In Fig. 2(a) we compare the exact scaling function for d = 3 (see Eq. (70)) 
to the Gaussian one (Eq. (49), formally d = A). In both cases X* =Q increases monotonically 
from the value 1 — 2/d for u — to d/(d + 2) for u — > oo. 

As in the case of the Gaussian approximation (see Sec. Ill), we consider the momentum- 
dependence of the FDR. As a function of the two scaling variables y = q z s = q 2 s and 
u = s/r m (where r m is given in Eq. (67)) one finds, from Eqs. (68) and (56): 



d-2 



+ y 



l l 



2 1 + u 



dx (1 — x) 



(d-4)/2 



U 



1 + u 



-x 



-2xy 



(74) 
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12 3 4 

y = q z s 

FIG. 3. Non-monotonic behavior of X™ as a function of y = q z s and different values of 2 < d < 4 
for the exactly solvable case n — ► oo of the O(n) model relaxing from a disordered initial state 
(mo = 0). The Gaussian result corresponds to d > 4. 

from which one recovers Eq. (70) for y = and the Gaussian result (50) for d — 4 and 
it — > oo, as expected. Taking into account that J^dx (1 — x)~ € ^ 2 e~ 2yx = 1/(2?/) + e/(8y 2 ) + 
e(2 + e)/(32|/ 3 ) + 0(y~ 4 ) one easily determines the behavior of for y >> 1: 

+ 0(y- 3 ). (75) 

Accordingly, for y = q 2 s ^> 1 the FDR converges to the asymptotic value 1, as heuristically 
expected on the basis of the fact that the mode actually responsible for the aging behavior 
is the homogeneous one q = 0. We note also that this asymptotic value is approached from 
above for < u < l/(y/2/e—l) and from below for u > l/(y/2/e—l). Given that X* < 1 for 
y — 0, X* has a non-monotonic behavior as a function of y for fixed < u < l/(y/2/e — 1) 
and in particular for u — 0, corresponding to an initially disordered configuration. This 
non-monotonicity disappears in the Gaussian limit d — > 4 and was already noticed at the 
first order in the e-expansion for the relaxation of the 0(n) model from a disordered state 
(see Fig. 1 in Ref. [7] and Ref. [3] for details). In fact, the expansion in e = 4 — d of Eq. (74) 
for u = agrees with the limit n — > oo of the results presented in Ref. [3] (cp. Eq. (3.14) 
therein). In Fig. 3 we provide a plot of Eq. (74), as a function of y — q 2 s for itlq = and 
d — 3, 3.5, and 4, the latter corresponding to the Gaussian model. The non-monotonic 
behavior is still quite pronounced for d — 3 and the overshoot decreases upon increasing the 
dimensionality d (cf. the dotted line corresponding to the unphysical dimension d = 3.5, 
reported for comparison), while shifting towards larger values of y. In the opposite limit, 
i.e., u — > oo (corresponding to r m = 0), a monotonic increase of X* as a function of y — q 2 s 
is expected. In Fig. 2(b) we compare the behavior predicted in this case by Eq. (74) (with 
u — > oo) for d — 3, to the Gaussian result (50) (formally, d — 4). We also indicate the 
limiting curve for d — > 2. In these cases X* increases monotonically from its initial value, 
given by d/(d+ 2), to the asymptotic value 1, which is approached according to Eq. (75). 
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V. ONE-LOOP FLUCTUATION-DISSIPATION RATIOS 



In the previous sections we have discussed the non-equilibrium behavior of models whose 
action is essentially quadratic in the fields (either because non-Gaussian terms were neglected 
or turned into effective Gaussian ones). In this section we account for the effects of non- 
Gaussian fluctuations within the perturbation theory up to one- loop (i.e., 0(e) in the e- 
expansion where e = 4 — d). In particular we compute the longitudinal and transverse 
correlation and response functions for a quench from a state with initial magnetization 
m — > oo to the critical point. As we have argued from general scaling arguments, a finite 
mo gives only corrections to the leading long-time behavior. The first step in this calculation 
is the solution of the one-loop equation of motion, presented in the next subsection. Then 
we calculate the connected (longitudinal and transverse) correlation and response functions 
(and the associated FDRs) in the following subsections, reporting the relevant details in 
Appendix B. 

For future reference we recall that the standard static exponents /3, 5, v of the universality 
class we are presently interested in are such that /3/u = (d—2+r))/2, 5 = (d+2—i])/(d—2+i]) 
with r] = 0(e 2 ) and z = 0(e 2 ) (see, e.g., [34]), and therefore [cp. Eqs. (38) and (39)] 

a = 0(e 2 ), e a = ^ = ~ + U0(e 2 ), and W = - A = _i + £ + (e 2 ) . (76) 

A. The equation of motion and its solution 

At one-loop level the tadpole contributions to the equation of motion (28) are 

I a (t) = J (dg)C£°(f , t) = 2N d r a t 1 - d l 2 , (77) 
Ut) = J (dq)C^°(t, t) q = 2N d r^- d / 2 , (78) 
where N d = 2/[(4ir) d / 2 T(d/2)}, and 

96 

Tu = (87r)^/ 2 (8 - d)(6 - d)(4 - d)(2 - d)N d ' (?9) 
4 

T " = (8n) d / 2 (4-d)(2-d)N d ' (80) 
Accordingly, the one-loop equation of motion is [see Eq. (28)] 

= d t m + ^m 3 + g N d (V CT + ^^r^j t^m = +O(g 2 ) (81) 

= dtm + I m 3 + g Qri £-*l* m + O(g 2 ) , (82) 
where g = N d g and 

n — 1 , 00 , 

rd = r a H — r„ , (83) 
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(2) 



1 



/ 

* / 
< * < 



(3) 



(4) 



FIG. 4. One-loop diagrams contributing to the longitudinal response function R^(t,s). The 
corresponding expressions (reported in Appendix B 1) are indicated by Ii(t, s), i = 1, . . . 4. 



which is solved by 



m(t) 



1 - 



9o 



3-d/2 



6/2 



+ 0(g 2 ) ■ 



(84) 



Expanding in e = 4 — d, one obtains 



m(t) = (l-g ^j 1 -g (^lnt + c - ^) + O(e 2 ,eg ,gl)j , 



(85) 



where c_i = — (n + 8)/ 12 and c = 9/16 — (n + 8) (In 2 + 7^)/24 are the coefficients of the 
expansion of = C-i/e + Co + 0(e). The dimensional pole ~ 1/e can be cancelled out by 
a proper renormalization of the bare magnetization (see, e.g., Ref. [34,32,33,30] for details). 
Once the renormalized m{t) is expressed in terms of g = g + 0(g 2 ), its scaling behavior 
at the RG fixed point g* = 6e/(n + 8) + 0(e 2 ) clearly shows up: m(t) ~ + 0(e 2 ) with 
q = 1/2(1 + c_ 1 g*) + 0(e 2 )=P/(isz). 



B. The Longitudinal response function 



At one-loop order the expression for the response function gets modified because of the 
one- loop term contributing to the magnetization m(t) in Eq. (85), and therefore 



= exp 



-q 2 (t-s) 



j\t'm 2 (t') 



S\3/2 

t. 



i ~ 3 t e/2 -s e/2 ^,_ 2 ; 



(86) 
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In addition, the interaction vertices of L\ [see Eq. (16)] yield four contributions which are 
depicted in Fig. 4. In terms of those diagrams the response function is given by 



Rl(t, s) = i2J°(f , s) - f/, - go^h + 9oh + 9o^h + 0{gl) . (87) 

The relevant zero- momentum integrals I iy % = 1, . . . , 4, have been calculated in Appendix B 1. 
Collecting the corresponding expressions we obtain, in terms of the renormalized coupling g 
and of the ratio x = s/t (0 < x < 1), 

- 2 61 + 2n 



R^ =0 (t,s)=x^ 2 
where 



l + g 



n + 8 
24 



In x + — 



7T 



24 



+ 0(g 2 } eg), 



26 + n 
12 



l + ln(l -x) 



- 1 



x 



-x 



;Li 2 (^) , 



(89) 



with /£(0) = 0. At the fixed point g* = 6/(n + 8)e + 0(e 2 ), Eq. 
expected scaling form (34) with 

- 2 61 + 2n" 



7T 
T 



24 



+ 0(e 2 ), 



F^x) = l + ~g*r R (x) + 0(e 2 ). 



can be cast in the 

(90) 
(91) 



C. Longitudinal Connected Correlation function 

As in the case of the longitudinal correlation function gets a one-loop correction 
Cq'°(t, s) due the one-loop term in the magnetization m(i) [see Eq. (85)]. This correction 
has been evaluated in Ref. [30], Eq. (B.29). In making use of that result we have to keep in 
mind that the factor in Ref. [30] accounts only for the contribution of longitudinal modes 
(i.e., rd would be r a in the current notation) whereas here it contains also the contribution 
of transverse modes [see Eq. (83)]. The additional interaction vertices contribute via the 
Feynman diagrams depicted in Fig. 5, in terms of which 

(%(t,8) =C?(t,8) + ^h + g ^I 8 +g I 9 +g ^I 9 -^ h~ go^ho + 0(g 2 ) . (92) 

Using the results reported in Appendix B 2 we obtain the zero-momentum longitudinal 
correlation function {x = s/t) 



C° =0 (t,s)= 1 -sx 3 / 2 



1 ~( ™ + 8 , 377r ' 25n 18913 \ _. a . , , _ 2 

1 + °(—2T lnx + W " 432 ~ W ) + gf ^ s/t) + °^ e9) 



(93) 

which has, at the RG fixed point, the expected scaling form with the correct exponents and 

* /37tt 2 25n 18913\ 9x , 

2A a r = l + o* + 0(e 2 ) , 94 

c y V 160 432 8640 J V ; ' K J 

FZ(x) = l + ~g*fc(x) + 0(e 2 ), (95) 
where fc(x) is given by Eq. (B31). 
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FIG. 5. One-loop diagrams contributing to the longitudinal correlation function Cq(t, s). The 
corresponding expressions are indicated by Ii(t, s), i = 5, . . . 10 and are computed in Appendix B 2 



D. Longitudinal FDR 



It is now easy to compute the longitudinal FDR from Eqs. (40), (38), (90), and (94) 



X? 



A£(l-0 ff ) 5 



4 1895 + 76n - 162tt 2 ^, 2N 

-e + 0(e 2 ) 



1800(8 + n) 



(96) 



which for n — 1 reproduces the result of Ref. [30] for the Ising model. Close to d — 4, X™ 
(n > 2) decreases as the dimensionality decreases. A monotonic behavior has to be expected 
down to the lower critical dimension d = 2, where it is known that X%° = 1/2 [24]. In the 
limit n — > oo the FDR becomes 



X~(n 



°° ) = J-4i) e+0(e2) ' 



(97) 
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FIG. 6. Left: Longitudinal FDR X£° as a function of the dimensionality d. The unconstrained 

([u], see Eq. (96)) and constrained ([c], see Eq. (101)) estimates for n = 2,3 and oo are shown 

together with the result X°°(sph) of Ref. [26] for the spherical model (dashed line). The dashed 

straight line indicates the expansion of X°°(sph) close to d = 4, which does not coincide with the 

limit n — > oo of Eq. (96). Right: Transverse FDR as a function of the dimensionality d. The 

unconstrained ([u], see Eq. (115)) and constrained ([c], see Eq. (118)) estimates for n = 2, 3 and oo 

are shown together with the exact result (71) for n = oo (dashed line). For comparison we show 

also the results of the e-expansion (e = d — 2) of Ref. [24], Eq. (73), for n = 3 and n = oo. 



which differs from the e-expansion of the result reported in Ref. 
for the spherical model: 



(see Eq. (8.116) therein) 



X°°(sph) = - 
5 



21n2 
15 



53 



7T 



900 120 



e + 0(e 2 ). 



(98) 



For finite n, a first rough numerical estimate of X%° in the interesting and non-trivial physical 
dimension d = 3 can be obtained by evaluating equation (96) for e = 1. For instance we 
have ([u] denotes unconstrained estimates, to distinguish them from the constrained ones 
introduced below) 



X?[u](n = 2) = 0.775, 
X?[u\(n = 3) = 0.774. 



(99) 
(100) 



However, the reliability of such results can be questioned, both because they are based on a 
one-loop approximation and mainly because for e = 2 one would find ^°[w] ~ 0.75, which 
is quite far from the exact result 1/2. To obtain more reliable estimates without computing 
the 0(e 2 ) contribution (which seems to be a difficult task, requiring the calculation of almost 
fifty Feynman diagrams) one can constrain [c] the O(e) result to assume the exact known 
value at d — 2, as usually done for this kind of expressions (see, e.g., [20,40]). Assuming a 
smooth behavior in e up to e = 2 one gets 



+ 



3 

10 



1895 + 76n - 162tt 2 3 
1800(8 + n) 20 



(101) 
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FIG. 7. One-loop diagrams contributing to the transverse response function K£(t,s). The 
corresponding expressions, indicated by ij(t, s), i = 11, . . . , 14, are calculated in Appendix B3. 



which has the same e-expansion as Eq. (96), assumes the exact value in d = 2 and there- 
fore it is expected to converge more rapidly to the correct three-dimensional result. The 
estimate (101) yields, for e = 1, 

X™[c}(n = 2) = 0.713, (102) 
X™[c](n = 3) = 0.712, (103) 

which are considerably smaller than the corresponding unconstrained estimates. A more 
robust field-theoretical estimate would require the knowledge of higher-order terms in the 
e-expansion. In Fig. 6 we compare the unconstrained [u] and constrained [c] estimates 
for n = 2, 3, oo as functions of the dimensionality d ranging between the lower and the 
upper critical dimensions of the model (2 and 4, respectively). For comparison we report 
the exact result X°°(sph) for the spherical model (dashed line), given by Eq. (8.113) in 
Ref. [26] (cp. also Fig. 4 therein). As already pointed out, our result (97) close to d — 4 (see 
the straight line n = oo[u] in Fig. 6) differs from the corresponding one (98) obtained in 
Ref. [26] (represented as a straight dashed line in Fig. 6). On the other hand, the numerical 
discrepancy between these two expressions is quite small, resulting in a constrained estimate 
X%°[c](n = oo) which differs from X°°(sph) by at most 3%. 



E. Transverse response function 



The transverse response function, like the longitudinal one, gets modified because of the 
one-loop contribution to m(t): 



<fo(M)=e-^- s H/>W) 



1 + ^0 



t e/2 _ s e/2 



l + e/2 e/2 



-q 2 (t-s) 



(104) 
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In addition, the interaction vertices in Ci [see Eq. (14)] yield the diagrammatic contributions 
depicted in Fig. 7, in terms of which the transverse response function reads 



Rl^ s) = RZ°(t,s) - f/ n - ^ 



(105) 



6 ^ 6 9 ° 9 

Using the expressions of ij, % = 11, ... , 14 reported in Appendix B3 we find for q = 

(x = s/t) 



T31T 

rt q=0 



(t,s) 



X 



1/2 



9 



where 



ln(l 



x) 



4a; 



— x ln(l — x) H x ln(l — x) 

24 3 V ; 12 v ; 



(106) 



(107) 



(/r(0) = 0). Equation (106) displays, at the RG fixed point g = g*, the correct scaling 
behavior [see Eq. (36)] with the proper exponents [see Eq. (76)] and 



An 



l -k +0{t2) 



FZ(x) = l + ffZ(x) + 0(e 2 ). 



(108) 
(109) 



F. Transverse correlation 



As in the case of Cq, the transverse correlation function gets a one-loop correction 
Cq'°(t,s) due the one-loop term in the magnetization m(t) [see Eq. (85)]. This correction 
can be computed from Eqs. (24) and (104), finding 

C^ =0 (t,s)=2 f S dt'Rl£ (t,t')Rl£ (s,t') 
Jo 



= sx 



,/2Ji +*w ,/2 



x- e / 2 - 1 2 
+ 



6/2 



2 + e/2 



O(~9 2 o) ■ 



(110) 



The additional vertices contribute via the Feynman diagrams depicted in Fig. 8, according 
to 

C-(t, s) = C^(t, s) - g ^I 15 - |/i6 + f h7 + |/is + f /i9 + 0{gl) . (Ill) 

Using the results reported in Appendix B4 for ij, i = 15, . . . , 19, we obtain the zero- 
momentum transverse correlation function 



C*(t, S ) q=0 = 8 (0 



S\l/2 



n + 8 , 127 n . . > % ., 

1 + /; ( ___ lnx + _ + _ + ) + {g\ eg) 



;il2) 



which displays the expected scaling behavior [cp. Eq. (37)] with the proper exponents [see 
Eq. (76)] and 



127 n 



AZ = l + g* [— + —) + 0(e 2 ) , 
c y \432 24/ v ; ' 

FZ(x) = l + ~g*fZ(x) + 0(e 2 ), 



(113) 
(114) 



where fc(x) is given by Eq. (B53). 
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FIG. 8. One-loop diagrams contributing to the transverse correlation function C™(t,s). The 
corresponding expressions, indicated by Ii(t, s), i = 15, . . . , 19, are reported in Appendix B 4. 



G. The transverse FDR 



The expressions of the transverse response and correlation functions allow the determi- 
nation of the FDR according to Eq. (40) [see also Eqs. (108), (113), and (39)] 

X - -^(l-^)-3"l08(8 + n) e + ° (e) ' (U5) 

which is a decreasing function of e, a feature that we expect to persist down to the lower 
critical dimension (d = 2 for the case of interest n > 2). Note that the coefficient of the 
correction of 0(e), i.e., 1/18 + l/[108(8 + n)], changes less than 2% with n > 2 and therefore 
X£°, at least close to d — 4, is practically independent of the actual number n — 1 of 
transverse components. In the limit n — > oo, reduces to 

X-(n = oo) = ^-l e + 0( e 2 ), (116) 
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in agreement with the e-expansion of the exact result X£°(n = oo) = djid + 2), obtained 
in Sec. IV. As for finite n, we proceed as in the case of the longitudinal FDR, first by 
giving unconstrained estimates (obtained by setting e to the proper value in Eq. (115) and 
neglecting 0(e 2 ) corrections) and then improving them by taking advantage of the known 
value for d = 2. The direct estimates for the physical dimension d = 3 are (they differ on 
the fourth digit) 

X^>](n = 2,3) =0.610. (117) 

These unconstrained [u] estimates are expected to work better than in the longitudinal case, 
because for e = 2 (i.e., d = 2) they give X^[u](n = 2, 3) ~ 0.55, which is only 10% above 
the exact value 1/2. To confirm this expectation we compute the constrained FDR 

A?H = i + (l-i) 

which gives, in three dimensions, 

X™[c](n = 2,3) = 0.597. (119) 

This estimate is indeed quite close to the unconstrained one, suggesting that in this case 
the one-loop expression provides more reliable results than it does for the longitudinal FDR. 
In Fig. 6 (right panel) a comparison between the various estimates as functions of the 
dimensionality d is provided for 2 < d < 4. The unconstrained ones for n = 2, 3, and oo 
(denoted by n — 2[u], n — 3[u], and n — oo[u], respectively) almost coincide due to the weak 
dependence on n of the expression (115). As a consequence, the corresponding constrained 
estimates (n = 2[c], n = 3[c], and n = oo[c], respectively) are not distinguishable on this 
scale. For comparison we report (dashed line) the exact prediction in the limit n — > oo 
(see Eq. (71)) and the prediction (73) for the e-expansion (e = d — 2) close to the lower 
critical dimension for n = 3 and n = oo (denoted by n = 3[e] and n = oo[e], respectively). 
As already pointed out in Sec. IV the result n = oo[e] is indeed tangential to the exact 
result n = oo, exact (dashed line, Eq. (71)). The constrained estimate n = 3[c] provided 
by Eq. (118) does not account completely for the behavior of X^{n = oo) close to d — 2, 
as predicted by Eq. (73), but reproduces only its value for d = 2. Reliable interpolation 
formulas accounting properly for both the behaviors close to d = 2 and d = 4 require the 
knowledge of higher order terms in e and e than those currently available. 

VI. CONCLUSIONS 

In this paper we have studied the non-equilibrium dynamics of a system belonging to the 
0(n) universality class and relaxing at criticality according to a purely dissipative dynamics 
from a state with initial magnetization M Q (induced via an external magnetic field). In the 
aging regime (t ^> s ^> r m ~ M 1 ^ K ) we found that the zero-momentum transverse (tt) and 
longitudinal (a) response and connected correlation functions have the scaling forms 



1 

6 



/ 49 + 6n 
V 108(8 + n) 



1 

12 



(118) 
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^=o(*> s ) s) a {t/s) e °F° R (s/t) , (120) 

q =0 (i, s) = ^ sit - snt/sf'Fgis/t) , (121) 

<=o(^ s) = (t - sfit/s^F^s/t) , (122) 

C£ =0 (i, s) = AZ s(t - sY{t/s)^F^{s/t) , (123) 

independently of the actual value of M ^ 0. We argued that the exponents 9 a ^ and 9'^ 
are not new non-equilibrium quantities, but they can be written in terms of equilibrium 
exponents as 

0, = £ = (124) 

vz 

= = (125) 

uz 

The two response and correlation functions define two (in principle different) FDRs and 

7T 

We solved exactly the dynamics of the transverse modes within the Gaussian approxima- 
tion and we derived all the relevant universal ratios and scaling functions (the dynamics of 
the longitudinal modes was already discussed in Ref. [30]). In particular we found that the 
FDRs of longitudinal and transverse modes are different, making the effective temperature 
(defined via X°°) observable-dependent already within the Gaussian approximation. Then 
we solved exactly the dynamics of the transverse modes in the limit n — > oo. The resulting 
FDR interpolates correctly (as a function of the aging time) from the well-known spherical 
model value X°° = 1 — 2/d [8] at M = and the new result 

= -J— , (126) 

at any M ^ 0. This result coincides with the one for local spin observables in the spherical 
model [26]. 

In order to provide more accurate predictions in physical dimensions d — 2, 3 for the 
physical models with n = 2, 3 we computed the first-order correction in the e-expansion to 
the equation of motion, to the response and correlation functions and from them we obtained 
the FDRs. In particular for the longitudinal modes we found 

4 _ 1895 + 7ft,- 162,' ^ 

5 1800(8 + n) K J K J 

which differs for n = oo from the e-expansion of the corresponding spherical model result [26]. 
To our knowledge this is the first example of a universal quantity taking different values in 
these two models, which are known to be equivalent in equilibrium [41]. It is not yet clear to 
us whether this discrepancy is actually due to an effectively different initial state or to non- 
Gaussian corrections to the action which have to be taken into account for the longitudinal 
mode (see Ref. [26]). For comparison, it would be interesting to calculate the leading non- 
Gaussian contributions to the 0(oo) model. For finite n we provided an estimate of by 
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constraining it to assume the known value in two dimensions (see Eq. (101)), resulting in 
X£°(n = 2, 3) ~ 0.71 in three dimensions. For the transverse FDR we obtained 

2 49 + 6n 2 
X - -3"l08(8 + n) e + ° (e) - (128) 

The constrained estimates (see Eq. (118)) in d = 3 yield X£°(n = 2, 3) ~ 0.60. 

An interesting extension of the work presented here would be to consider the non- 
equilibrium evolution at criticality of the 0(n) model starting from an equilibrium magne- 
tized state which has been prepared in the low-temperature phase T < T c in zero magnetic 
field. In this case the transverse modes are critical even in the initial state, affecting in a 
non-trivial way the ensuing evolution. 
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APPENDIX A: DETERMINATION OF THE h FUNCTION 

In this appendix we determine the function h(x) [introduced in Eq. (60)] which charac- 
terizes completely the dynamics of the transverse modes in the limit n — > oo. 

First of all we express the self-consistency condition Eq. (57) in terms of h only. To this 
end we note that Eq. (60) implies: 



and that, from Eqs. (62) and (56), 



/ ^c^^jy^gi j (d9) e-» <f - f -,. 2iVdA ^' dt ^ Aj(2A2(i _ i0) , 

(A2) 



where we have introduced the kernel 



2^ /A2 



K d {r) = N-^- d J (dg)e~^/ A . (A3) 

Its explicit expression depends on the specific implementation of the large-momentum cut 
off. For later convenience we have introduced the factor N d = 2/[(47r) d//2 r(<i/2)] in the 
definition of K d . Then we compute 

/ A (dg)C-(oo, oo) = [\dq)^-— = N^I^/ \ 2 ) , (A4) 
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where I d {z) is given by 

I d {z) = N-'A-^ fidq)^!^ (A5) 

(I d (0) is finite for d > 2). Using these results and Eq. (61), the consistency condition (57) 
can be expressed in terms of h only 

A W) = r - + -6~ A \J dT W) ' d{ " } " d( 007 } J + T fcOO ' m 

where r = 2A 2 t and r' = 2A 2 t' are dimensionless time variables. This equation turns into a 
linear differential equation for h: 

h'{r) = ^h(r) + u y^dr' h{r')K d {r - t>) - 7 d (r 0O /A 2 )/ i (r) } + , (A7) 

where the dimensionless coupling constant is given by u = gN d A d ~ 4 /Q and which can be 
easily solved by using the Laplace transform g(s) = C[g{r)] = J °°dr(y'(r)e~ s ' r , defined in the 
half-plane Re(s) > 7, where 7 depends on the analytical properties of g. Once g(s) is known, 
the original function is obtained via the inverse transform g(r) = (27ri) _1 f^,**™ ds g(s)e ST 

for 7' > 7. Transforming Eq. (A7) one gets C[h'(r)] = —h(0) + sh(s) where, from the 
definition (60), h(0) = 1/2, and 

/ dre" ST / dr'h(T')K d (r-T') = / dr'e- ST ' h(r') dre' 8 ^-^ K d (r - r') = h(s)K d (s) , 
./o ./o ./0 Jt' 

(A8) 

where (see Eq. (A3)) 



dre—^(r) = A"" jf dre— J {dq)^'^ = A~ d J (dg) q2/A2 + g = W) 

(A9) 

Accordingly, Eq. (A7) becomes 

I ^ V *• 1 777 ^ 1 

"2 + SMS) = ^ h{s) + " [/d(s) ~ /rf(rWA2) ] ^ + 6 ' (A10) 
whose solution is 

1 lm 2 1 

h(s) = ; 1 6 A ' S . (All) 

* - ^ - - ^oo/A 2 )] 

First we note that, as expected, h(s) has a pole for 3 = 300 = f^/A 2 independently of the 
value of u. For u — 0, /i(s) is analytic in the complex plane except for the poles at s = 
and s = 0. Accordingly, the corresponding inverse Laplace transform is 

1 1 m 2 1 

fc(r) = -e'~ T + -^— (e'~ T -l), (A12) 
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7' 
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FIG. 9. Analytic structure of h(s) and contour of integration for the computation of its inverse 
Laplace transform. Poles are indicated by a cross, whereas the branch cut by the thick solid line. 
7 is the original contour of integration to compute the inverse transform of h(s), whereas 7' is its 
convenient deformation in the case Sqq = (see Appendix A 1). 



i.e., 

h(2A 2 t) = + (e 2 ^* - 1) , (A13) 

which is actually independent of A, as expected. From this expression one recovers all the 
Gaussian results of Sec. Ill in the limit — > (corresponding to the critical point). 

Because of Id{s), h(s) has a branch cut on the real axis for s_ < Re(s) < and u > 0, in 
addition to the poles for s = and s = 0. If s_ is finite then an isolated pole s_oo is also 
present on the real axis with Re(s_oo) < s_ (see Appendix Al). In the scaling limit we are 
interested in, A 2 t > 1 (r > 1), the relevant contribution in the inverse Laplace transform 
comes from the singularities with the largest real part. If > 0, this is given by s = Sqo, 
leading, as in the Gaussian case, to an expected exponential decay of h~ l {r). For foo = 0, 
instead, 

1 1 ml 1 

h( s ) = \ I , n , (A14) 

and the leading singularity comes from the branch cut. Using the asymptotic expansion for 
Id(s) derived in Appendix A 2 we can write the denominator of Eq. (A14) as 

00 

s - u[I d (s) - I d (0)} = a Q s d ' 2 ' 1 + 6 s + Ws d/2 ~ l+k + b k s 1+k ] , (A15) 

k=i 

where 

a Q = uK d and b = 1 — uctd, (A16) 

and Kd = 7r/[2 sin(7r(d — 2)/2)] is a universal coefficient, i.e., independent of the actual 
regularization procedure, whereas ad does depend on it (with the universal feature a4- e ~ 
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1/e for e — > 0). Accordingly, for s — > and 2 < d < 4 (so that e = 4 — d > 0) one has 

1 



s-u[I d (s)-I d (0)] 



= a s 



l„-l+e/2 



1 



E 

m,n=0 



m+ne/2 



(A17) 



where Co 5 o = 0, Co jra = (— &o/ a o) n an d & h the other coefficients can be computed in terms of 
at and bk- [Note that for d > 4 the leading behavior for s — > is, apart from a coefficient, 
the same as in the case u — 0, leading to the expected Gaussian-like behavior.] The inverse 
Laplace transform of h(s) for r ^> 1 is dominated by the behavior of h(s) for s — > 0, and it 
can be computed following a suitable contour of integration in the complex plane (indicated 
as 7' in Fig. 9) on the two sides of the branch cut. Recalling that £ _1 [l/s a ] = T a ~ 1 /T(a) : 



one concludes that 



Kr»l) = \ 



.-6/2 



LfJ-i 



r(l-e/2) 



CQ,n 



-(n+l)/2e 



n=l 



(l-(„ + l)/2c) 



(A18) 



+ 



a m, 



o 



6 A 2 
+ 0(r- 1 ). 



.l-e/2 



LlJ-i 



+ E 



l-(n+l)/2e 



LlJ-i 



r(2 - e/2) ^ r(2 - (n + l)/2e) ^ T(l - (n + l)/2c 



-(n+l)/2e 



The confluent corrections ~ r~ ne l 2 (n > 2) to the leading scaling behavior vanish identically 
if c 0i n = 0, i.e., for b = 0. According to Eq. (A16) this occurs if u — u* — a^ 1 , which is 
indeed the fixed-point value for u (note that if dimensional regularization is employed, this 
leads to the well-known result u* — e). 

Accordingly, the scaling limit (i.e., A — > oo with rn^t fixed) for 2 < d < 4 is characterized 
by the fixed-point expression 



h(2AH > 1) 



(2A 2 t) 



-6/2 



u* 2r((d-2)/2) 
which is the result reported in Eq. (63) 



1 + 



3(d-2) 



rn^t + 0((2A 2 t)- £ / 2 ) 



(A19) 



1. Poles of /i(s) 

Here we investigate the properties of analyticity for s e C of 

1 lmgl 

^) = — ~ 2+ ^7\~ s T ( s\ ' ( A2 °) 

s - p - u[I d (s) - I d (p)\ 

where p = foo/A 2 > 0. Clearly, the numerator has a pole for s = whereas the denominator 
vanishes for s = p, independently of u > 0. id(s) [see Eq. (A5)] has a branch cut on the 
real axis for s_ = — q„ iax /A 2 < Re s < 0, where q max is the maximum momentum in the 
integration domain (|q max | ~ A). Depending on the actual definition of the momentum 
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integral J (dq), s_ might be not finite. Possible additional poles are given by the solution 
of 

s-p-u[I d (s)-I d (p)]=0. (A21) 
Note that, as a consequence of Eq. (A5), 

W - Up) = - rt^-A^/w w + A , t) 1 w + AV) • (A22) 

and therefore additional poles are given by the solutions of 

J<( ,),^/;) (q , +ft) ' tf + Ar -l. (A23) 

It is easy to see that the imaginary part of the integral vanishes only on the real axis 
and therefore possible poles have to be located there and we can restrict our search to 
s G (-00, s_) given that J d (s G R + ) > 0. Note that J d (s < s_) < 0, J d (s — > -00) = 0", and 
that J' d {s) < (as it is clear by differentiating the integral). In addition, J d (s — > s_) diverges 
logarithmically and therefore Eq. (A23) admits a unique solution s_oo < s_ for 11 ^ 0. In the 
case of a regularization with a sharp cut-off A one finds s_ = — 1 and it is easy to determine 
the asymptotic behavior of J d (s) for s — > — 1~: J d (s — > — 1~) = — (l+p) -1 In (— 1 — s) + finite 
and therefore Eq. (A23) admits a unique solution s_oo < —1 for « ^ 0. In particular, for u 
small one has s„„o ~ — 1 — q~ 2 ^ 1+ p)/ u . 



2. Expansion of I d (s) 

Let us consider a possible definition of the regularized integral, via a cut-off function 
(p(x) (defined for x e i.e., 

/ A (dg) ^ / (dg) v.(q 2 /A 2 ) , (A24) 

where we assume that tp is a smooth function which does not change the small-momentum 
behavior of the integrand, i.e., tp(0) = 1 and decays exponentially fast for x — > 00. For 
later convenience we also assume that <p(x) can be expanded around x = (or, at least that 
(f'(0) is finite). These properties are also encoded in the analytic properties of the Mellin 
transform <p(s) of (fi(x), defined as 

POO 

<p(z) = / dxx'-^x) , (A25) 







which converges in the complex plane for Re(z) > (in case of an algebraic decay of (p(x) 
for x — > 00 this domain becomes a strip). Its analytic extension to z G C is characterized by 
poles on the real axis for the non-positive integers n = 0, -1, -2, . . . with residues <^ (n) (0)/n!, 
due to the small-x behavior of <p(x). The original function (p(x) is recovered by computing 



1 



c+ioo 



V(x) = — dz<p(z)x- g , (A26) 
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with c > 0. 

Using this representation of tp in the definition of I d [see Eq. (A5)] we find 

1 pc+ico 

U s ) = — dz^s^ 2 - 1 -* f- -, (A27) 

dK } 2mJ c _ lQO 2sm^^(d/2-z) , 1 ; 

where d/2 — 1 < c < d/2. This expression allows us to determine the asymptotic behavior 
for s — > 0, by shifting the integration contour towards negative values of Re z: 

I d (s) = -<p{d/2 - 1) + 2sin( 7r 7rrf/2) ^(0)^ /2 ' 1 - 2^(^/2 -2) S + E [A k s d ^ k + B k 8 1+k ] . 

(A28) 

Note that the first term, I d (0) (d > 2), depends on the specific choice of the cut-off function 
<p(x), whereas the second term is actually universal provided that y?(0) = 1. Higher-order 
terms do depend on ip(x). Let us consider the linear term is s. The analytic extension of 
<p(z) to the region — 1 < Re(z) < is readily computed as 



and therefore 



POD 

(z) = / dxx^lipix) -ip(0)}, (A29) 
Jo 



1 1 f°° 

a d = --^{d/2-2) = --J^ dxx d / 2 - 3 [ V (x)-l] (A30) 

behaves as ^ 

a 4 _ e = — h finite for e — > . (A31) 

Accordingly, Eq. (A28) can be cast in the form 

oo 

m = W) - «d s d/ ^ + a d s + J2 [A k s d ' 2 - 1+k + B k s 1+k ] , (A32) 



k=i 



where n d = 7r/[2 sin(7r(rf — 2)/2)], which leads to Eq. (A15). The same conclusion can be 
drawn by adopting different regularization schemes (see, e.g., Appendix B of Ref. [42]). 



APPENDIX B: ONE-LOOP FEYNMAN DIAGRAMS 



As usual in the computations of Feynman diagrams, it is convenient to calculate first the 
corresponding one-particle irreducible (1PI) parts. For the response functions they are given 
by the two tadpoles I a , 1^ (see Eqs. (77) and (78)) and by the four "bubbles" which can be 
written as (hereafter the response and correlation functions appearing in the integrals are 
the Gaussian ones) 
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(dg)i^(t, s)C*(i, s) = jf dxx 3 ^ - x)- d ^ = j^H L (t, s) , (Bl) 

B TL (t, s) = j (dq)Rl(t, s)C°(t, s) = - s B LL (t, s) = ^-J-H^t, s) , (B2) 

(dq)R^(t, s)CZ(t, s) = J dxx(t - x)- d ' 2 = J^H T (t, s) , (B3) 

B LT (t, s) = J (dq)Rl(t, 8 )C%(t, s) = S -B TT (t, s) = -J^ S -H T (t, s) , (B4) 



-d/2 



where 

H L (t,s)=r 3 [ dxx 3 (t-xY 
Jo 

= A L r 1+ ^ 2 - A L r\t - sy' 2 + B L r 2 s(t - s f 2 + c L r 3 s 2 (t - sf' 2 

+D L r 3 s 3 (t - s)- l+e/2 , (B5) 

H T (t,s)=t- 1 [ dxx{t-x)- d ' 2 = 
Jo 

= A T t~ 1+e ' 2 - A T r\t - sf/ 2 + B T t-\s(t - s)- l+e ' 2 , (B6) 

with 

6 

AL = (2 + e/2)(l + e/2)e/2(-l + e/2) ' (B7) 

6 

Bh = "(2 + e/2)(l + e/2)(-l + e/2) ' (B8) 
3 

° L = "(2 + f /2)(-l+e/2) ' (B9) 

^ = -rnwy (B10) 

^ = t /2(-l +t /2) - (BU) 

5 T = — — . (B12) 

(-l + e/2) V ; 

This parameterization of the 1PI parts of the integrals is particularly convenient, as explained 
in Ref. [30]. 

In what follows we shall make extensive use of the results of Ref. [30]. However, we alert 
the reader to the fact that the factor appearing in Ref. [30] has to be identified, in the 
current notation, with r a and not with = r a + (n — l)r„/3 (unless differently stated). 

1. Longitudinal response 

The integral l\ and 1% have been computed in Ref. [30], Eq. (A. 9) and (A. 11), respectively. 
J 3 is easily related to 1^. 

I 3 (t,s) = [ dt'R°{t,t')I„{t')R°{t',s) = —h(t,s) , (B13) 

J s r u 
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whereas J 4 is given by 



h(t,s)= I d*i / ds 1 R a (t,t 1 )V2m(t 1 )B T T(ti,s 1 )V2m(s 1 )R a (s 1 ,s) 

J s J s 

Q /s\3/2 r* /•*! 

/ dMi / ds!S^ 2 i7 T (ti,si) . 



(B14) 



According to the decomposition (B6) of Ht, this expression can be written as I4 = X^=i(°*)> 
where (a*) are easily calculable: 



(ai) 
(a 2 ) 



d*x / d Sl sftf = 1 ^ , g /oN - ^ 2 - S ^ 2 ) 



s(l + e/2) 

jf dt 1 J 1 ds lS ^(t 1 - Sl y^ 



l + e/2 



t 7T 2 



- - 1 - In - + 1 ( ( - - 1 j ln(t - s) - (1 + hit) In - + — - Li 2 (s/f) 



+ 0(e), 



Bt 



= f d*i dsi s^ti - Sl )- 1+£ / 2 = - In - + In tin - + U 2 (s/t) - ^ + 0(e) . 
Js J s ess o 



and whose sum is (x = s/t) 



6W 2 
(87r) d / 2 



X 



3/2 



--lnx+ 2 



In x ( 1 \ , . . 

I ( --1 )ln(l-x)+0(e) 



a; / 



(B15) 



2. Longitudinal correlation 

The integrals J 6 and I 7 have been calculated in Ref. [30], Eq (B.25) [see also Eq. (B.26)] 
and Eq. (B.27), respectively, whereas I 10 is trivially related to If. 



J-io — — h 



(B16) 



I g can be written as 

H rt 



I g (t,s)= j dt 2 j dt 1 R^ =0 (t,t 1 )V2m(t 1 )B TT (t 1 ,t 2 )V2m(t 2 )C^ =0 (t 2 ,s) + 
Jo Jt 2 

+ [ dt 2 [ dt 1 C^ =0 (t,t 1 )v / 2m(t 1 )S TT (t 2 ,t 1 )v / 2m(t 2 )^ =0 (s,t 2 ) (B17) 
Jo Jt 2 

which, using Eq. (B3), becomes 



I 9 {t, s) 



(st) 



-3/2 



s r t 
2 



dM 2 / dhhHxfa^) 
Jti 
t 



(87r) d / 2 

ft/ fT f S f S 

+s 4 / dt 2 t 2 - 2 / dti*if/r(*i,t 2 )+ / dM 2 / dhhHT^h) 

Js Jt 2 JO Jti 



(B18) 
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We denote the three terms in square brackets by (A), (B), and (C), i.e., [...] = [(A) + (C) + 
(C)]. In turn each of them can be written as the sum of the three contributions coming from 
the decomposition of H T according to Eq. (B6), leading to the following expression for I 9 : 



(87r) d / 2 



(B19) 



beA,B,C i=l 



The terms with % — 1 and 2 are the same as the corresponding ones calculated in Ref. [30], 
Eqs. (B.10), (B.ll), (B.15), (B.16), (B.20), and (B.21), with the proviso that A d therein has 
to be replaced by A T here (see Eq. (Bll)). We still need 



Bt 



B T s 4 
B, 



f ' dt 2 t\ f dhih-h)' 1 ^ 2 
Jo Jt 2 

[ dht^ 1 I dt 1 (t 1 -t 2 )- 1+e / 2 = -ln- + lntln- + Li 2 (s/t)-^ r + 0(e), (B21) 
J s Jto ess o 



(B20) 



PS l>S 

/ dhtf / d* 2 (* 2 -tl)- 
JO Jt x 



l+e/2 = S A 



>T JO Jt x 

Summing these nine terms we have 

3s 4+ ^ 2 (st)- 3 / 2 



1 25 1, , 

h - Ins +0(e) 

2e 48 4 ' 



(B22) 



where 



C(x) 



/9(t ' s) " (8tt)^ 
41n(l-x) 1 



1 - 2\nx \n 2 x 107 „, ' 



X 



3x 4 



ln(l - x) + x + — + — 



-ln(l-ar) 



(B23) 



(B24) 



[C(0)=0]. 

I 5 has been calculated in Ref. [30], Eqs. (B.5) and (B.6). In order compute I 8 we first 
need its 1PI part [we define t < = min{t, £'}] 

B T c T c {t,t>) = J (dq)[C-(t,t')] 2 = J (dqnttr'e' 2 ^ jf* dt ltl e 2 ^ jf* dt 2 t 2 e 2 ^ = 



4 (ft 
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dti I dt 2 
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1<<2 



(t + t'-h- t 2 y/ 2 • 



(B25) 



so that 
I* 



dn / dr 2 R^s,T 1 )V2m(r 1 )B^(n,r 2 )V2m(T 2 )R^t,r 2 ) 

PS f>T\ PS pt 

2 dri / dT 2 TiT 2 Bl^(Ti,T 2 ) + / dr x / dr 2 n^B^in, r 2 ) 

. Jo Jo Jo Js 



3(st)- 3 / 2 



3(st)-V 2 [2B + A] 



(B26) 
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Both A and B do not diverge for d — > 4. Therefore no dimensional regularization is required 
and one can calculate them directly in d — 4: 



/■S /■< />Tl /-Tl 

A = 4(8vr)^/ 2 / dri / dr 2 / dfi / dt 2 

JO </s JO JO 



(ri + r 2 - f i - t 2 



+ 0(e) 



= 4(87r)- d /V 



V(s/t) + hn2-^ 



+ 0(e), 
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PS PTl PT2 PT2 
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Jo Jo Jo Jo {n + r 2 - h - t 2 y 
7 1 



+ 0(e) 



4(87r)- d/ V ^--ln2 + 0(e), 
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where 



N 1 1 1 

^) = o + 7-T + 



1 W1 , / 1 5 1 1 1 

L W 1 — x) I I I 

8 ' 4x 3 ' Ylx 2 Ylx \ 24 24x 4 6x 3 Ax 2 Qx 



+ ln(l + x) -— - 



1 1 



1 ~r\ 



24 24x 4 6x 3 Ax 2 6x J 



[P(0) = 0]. Accordingly, 
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(?) 



3/2 



(B27) 



(B28) 



(B29) 



(B30) 



Summing all these contributions according to Eq. (92) one finds that the function intro- 
duced in Eq. (93) is given by 



Kit) = \a{x) + \b{x) + ^C(x) + ^V(x) , 



(B31) 



where A(x), B(x) are given by Eqs. (B.5) and (B.26), of Ref. [30], whereas C(x) and T>(x) 
are given by Eqs. (B24) and (B29) of this appendix. 



3. Transverse response 

The integrals for I n and I 12 are (see Eq. (A. 9) in Ref. [30]): 

I n (t,s) = / dt'R"(t,t')I a (t')R"(t',s) = -h(t,s) , 

J s S 

h 2 (t,s) = [ MR"(t,lf)I v (t')R ir tf,s) = —I 11 (t,s). 

J S r cr 

Ji3 is, instead, given by 

In(t,s)= I d*i / ds 1 i? 7r (t,t 1 ) v / 2m(t 1 )fi TL (t 1 ,s 1 )v / 2m(s 1 )i? 7r (s 1 ,s) 

J s J s 

6 /s\ 1 / 2 f f [ ll 

dMi / dsi Si H L (ti, si) , 



(B32) 
(B33) 



(87r) d / 2 \t 



(B34) 
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which is just I 2 (t, s)t/s, where I 2 has been calculated in Eq. (A. 11) of Ref. [30]. For J 14 we 
have 

Iu(t, s) = dii J d Sl IF(t, t 1 )v / 2m(t 1 )fi LT (t 1 , si)v / 2m(s 1 )i? 7r (s 1 , s) 



(87r) d / 2 \t 



s \l/2 /•* /•*! 

■ dMf 1 / dsiff r (*i,si), 



(B35) 



which can be split in three parts, according to the parameterization (B6) of Ht- The final 
result is 



J 14 (t,s) = — ^ (-) 



(87r) d / 2 Vi/ 



In x + In x + 



In X 7T 2 



+ 3(1 -x)-(l-x) ln(l - x) + 2Li 2 (x) + 0(e)] . 



(B36) 



4. Transverse correlation 



As usual, the tadpole integrals J15 and I±q are easy 

J 16 (t,s)= / dt' R"{t',s)I tT (t')C ir {t,t')+ I dt' R"(t\t)I a (t')C"(t\s) 
Jo Jo 



= 2N d r, S (i) 



1/2 



e/2 



2 + e/2 e 

^15(*,S) = —Iw(t,s) . 

There are two diagrams involving Blt(Ji, s\): 

h 7 (t,s)= [ dt 2 [ dt 1 R"(t,t 1 )V2m(t 1 )B LT (t 1 ,t 2 )V2m(t 2 )C"(t 2 ,s) 
Jo Jt 2 

+ [ dt! [ dt 2 C"(t,t 1 )V2m(t 1 )B LT (t 2 ,t 1 )V2m(t 2 )R 7r (s,t 2 ), 

Jo Jt x 

which can be written, via Eq. (B4), as 
6(st)" 1/2 



I 17 (t, s) 



(87T) d / 2 



(B37) 
(B38) 



(B39) 



dt 2 t\ [ dt 1 ti 1 H T (t 1 ,t 2 ) +s 2 [ dt 2 [ dht^Hrit^h) 

Jt 2 Js Jt 2 

+ / dht\ / dht^Hrit^ti) . (B40) 
Jo Jti 



In turn, we decompose each of the three terms in square brackets in three parts, according 
to the parameterization (B6) of H T . We do not report the details of the calculation, but 
only the final result 



I 17 (t,s) 



Qs 1+e/2 /s\V2 



(8%y/ 2 \t 



2 2, ln 2 x , 23 tt 2 , ,' 

lnx + -— + lnx + — - — + G{x) + 0(e) 

6 6 2 18 3 



, (B41) 
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where 



1 20 21n(l -x) +x 2 , M . nT . . . 
0(x) = ---—x- - — — + -xlog(l -x) + 2U 2 (x) 



3 9 

[0(0) = 0]. 

The expression for n 8 is 



ar 



(B42) 



hs(t,s)= [ dt 2 [ di lJ R 7r (t,t 1 )v / 2m(t 1 ) J B TL (t 1 ,i 2 )v / 2m(t 2 )C 7r (t 2 ,s) 
</o Jt 2 

+ / dti / dt 2 C 7r (t,ti)v / 2m(t 1 )5 TL (t 2 ,t 1 ) v / 2m(t 2 ) J R 7r (s,t 2 ), (B43) 
which can be written, using Eq. (B2), as 



hs(t, s) = 



6{st) 



(8%) d / 2 



[ dt 2 t 2 2 [ dtrt(ii,< 2 )+« 2 / dt 2 / dht^ 1 H L (h,t 2 ) 

JO Jt 2 Js Jt 2 

+ f dhtl [ dhq 1 H L {t 2 M) ■ 

Jo Jt\ 



(B44) 



We decompose again each of the three terms in square brackets in three parts, according 
to the parameterization (B5) of Hl. Also in this case we do not report the details of the 
calculations, but only the final result 



hs(t, s) = 



6s 1+e/2 /s\V2 



(8n) d / 2 \t 



2 2 



In x 3 



7T 2 17 



___ l n;r + __ + _ l nx + _____ + JF( X ) + 0(e) 



where 



x 2 lnfl — x) 
F{x) = 5 + - -41n(l -x) - — [ln(l - rr) + x] +6— - 2Li 2 (x) 



(B45) 



(B46) 



F(0) = 0]. 

There is a last 1PI diagram 



B T c L c (t,t') = J(dq)CZ(t,t')C°(t,t f ) 

\dg)4(tt')" 2 e- 2<?2(i+ *' ) dht\e 2 ^ [** dt 2 t 2 e 2 ^ 2 

Jo Jo 



4(W) 



/\-2 /•*< 



(8vr)^ 7 

in terms of which we express Jig as 



t< 



dt t / dt 2 



t?t 2 



(t + V-U-h)*/ 2 ' 



(B47) 



'■19 



/" dn f dr 2J R 7r (s,r 1 )v / 2m(r 1 ) J B^(r 1 ,r 2 )v / 2m(r 2 ) J R 7r (t,r 2 ) 
o Jo 



= 3(st) 



-1/2 



2 / dn / dr 2 S^(ri,r 2 )+ / dn / dr 2 B^(r u r 2 ) 

. Jo Jo Jo Js 



3(st)-^ 2 [2B + A] 



(B48) 
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A and B do not diverge for d — > 4 and therefore one can calculate them directly in d — 4: 



A = 4(8vr)- d / 2 / dn / dr 2 - / dt, / dt 2 

Jo is \ T 1 T V Jo Jo 



B 



t\t 2 



(ri + t 2 - ti - t 2 ) : 



+ 0(6) 



= 4(8tt) 



-d/2 2 



7T 



+ 0(6), 



(B49) 



PS r-T\ i />T2 l>T2 

= A{M)- d ' 2 / dn / dr 2 - / df! / dt 2 

Jo Jo (nr 2 ) J Jo Jo 



t\t 2 



(ti + r 2 - *i - t 2 y 



+ 0(6) 



4(8vr) 



-d/2 2 



29 7T 2 4 
48 24 ~ 3 



ln2 +0(e) , 



(B50) 



with 



£(x) 



dy 1 {y^x - 1) / dxi / dx 2 



x\x 2 



(1 + yi -xi-x 2 y 



and 5(0) = 0. Thus 



/iq = 



12 /s\V2 
-s 



24 



(B51) 



(B52) 



" 19 (8n) d /^\t- 

Summing all these contributions according to Eq. (Ill) (see also Eq. (110)) one finds that 
the function fg introduced in Eq. (112) is given by 



g(x)+F(x) + 2£(x) 
12 



(B53) 



where Q(x), J-{x), and £(x) are given by Eqs. (B42), (B46), and (B51), respectively 
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